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ABSTRACT 

Analyses have found a "haze" of anomalous microwave emission surrounding the 
Galactic Center in the WMAP sky maps. A recent study using Fermi data detected a 
similar haze in the 7-ray. Several studies have modeled these hazes as radiation from 
the leptonic byproducts of dark matter annihilations, and arguably no convincing 
astrophysical alternative has been suggested. We discuss the characteristics of astro- 
physical cosmic ray sources that could potentially explain this microwave and 7-ray 
emission. The most promising astrophysical scenarios involve cosmic ray sources that 
are clustered such that many fall within ~ 1 kpc of the Galactic Center. For example, 
we show that several hundred Galactic Center supernovae in the last million years plus 
a diffusion-hardened electron spectrum may be consistent with present constraints on 
this emission. Alternatively, it could be due to a burst of activity probably associated 
with Sagittarius A* occurring ~ 1 Myr ago and producing > 10 51 erg in cosmic ray 
electrons. Different models predict contrasting trends for the spectral index of the 
microwave and 7-ray spectrum as a function of angle from the Galactic Center that 
we show should be robust to cosmic ray propagation uncertainties. In particular, if 
the haze is from dark matter annihilations, it should have a very hard microwave and 
7-ray spectrum for which the spectral shape does not change significantly with angle, 
which we argue would be difficult to achieve with any astrophysical mechanism. Ob- 
servations with the Planck and Fermi satellites can distinguish between viable haze 
models using these signatures. 
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1 INTRODUCTION 



Using data f rom th e Wilkinson Anisotropy Probe (WMAP), 
iFinkbeinerl l|2004al ) discovered an approximately circular 
excess of diffuse emission centered around the Galactic 
Center of radial extent ~ 20° and luminosity Lh azc ~ 
10 37 ergs -1 between 20 and 60 GHz. Nicknamed the 
"WM AP haze," t his ra diation was the residual emission 
after IFinkbeinerl |2004a) subtracted templates for Galac- 
tic Ha, synchrot ron, and du s t emis sion from the WMAP 
foreground maps. IFinkbeinerl (|2004ah argued that the spec- 
trum of the haze was too hard to be due to synchrotron 
from cosmic ray electrons produced in supernovae. There- 
fore, this emission was i nitially attributed t o free-free emis- 
sion from > 10 5 K gas |Finkbeine r 2004a). However, in a 
follow-up study IFinkbeinerl (2004b) argued that the lack 
of diffuse X-ray emission from the Galactic Center ruled 
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out this hypothesis. If the haze is not due to free-free 
emissio n, it must be synchrotron. Interestingly, IFinkbeinerl 
(2004b) showed that the energetics and profile of the haze 
emission could be explained by synchrotron emission from 
the byproducts of dark matter annihilations for a generic 
weakly interacting ~ 100 GeV particle with a standard 
dark matter halo profile. This id ea has been pursued fur- 
ther in several sub s equent studies dHooper et al.ll2007l.[200; 



Zhan g et al.l 1200a ICholis et al.l l2009bl ; ICumberbatch et al 
20091 ; iHarding fe Abazauanll2010l ) ~ 

However, cosmic rays from astrophysical sources may 
also be able to create the haze emission. The total energy in 
synchrotron electrons required to produce the haze is 



^hazc 1Q 5i 



10 /iG 

B 



50 GHz\ 5 



erg. 
(1) 

Equation JT]) is derived by equating Lhaze with the syn- 
chrotron loss rate for an electron that emits at v in a mag- 
netic field of strength B, noting the relationship between 
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electron energy E and the synchrotron peak frequency v: 



E « 30 



50 GHz 



1/2 



10 /iG 



-1/2 



GeV, 



(2) 



where B± is the magnetic field component orthogonal to the 
particle motion. For comparison, ~ 10 erg is roughly equal 
to the kinetic energy of a single supernova explosion. This 
energy must be output within the time it takes an electron 
to radiate its energy via synchrotron and inverse Compton 
or 



At = 1.0 



B ± 50 GHz 



10 fj,G 



1/2 



Ur+U B 

10 eV cm _; 



Myr, 



(3) 

where U t and Ub are respectively the energy densities in 
radiation and magnetic fields. 

While it is conceivable that astrophysical sources could 
satisfy these energetics, the difficulties with creating the 
haze with astrophysical sources are threefold: 

(i) iDobler fc Finkbeinerl (|2008h find the haze spectrum 
in the microwave to be hard, measuring a power-law elec- 
tron spectrum with index 7 ~ 2. Such a spectrum is much 
harder than the galactic synchrotron towards other direc- 
tions (where 7 w 3). 

(ii) The distribution of known galactic cosmic ray sources 
appears to not be sufficiently concen trated toward the 
Galactic Center to creat e the haze |Finkbeinerl l2004bl ; 
IDobler fc Finkbeinerl 120081 ; IZhang et alj|2009h . The distri- 
bution of supernovae and pulsars is estimated to peak at 
« 4 kpc from the Galacti c Center, correspond ing to an an- 
gular separation of ~ 30° |Lorimer et al.ll2006t ). In contrast, 
the haze intensity is strongly increasin g with decreasing an- 
gle at < 20° from the Galactic Center ( Dobler fc Finkbeinej 
l200Sh . 

(iii) By construction, a similar excess towards the Galac- 
tic Center is not present at 408 MHz. If the haze is from 
synchrotron, this suggests that the source of the haze does 
not contribute significantly to the population of cosmic ray 
electrons at ~ 5 GeV. In addition, the sky at 408 MHz is 
dominated by emission from the galactic disk, without a sig- 
nificant enhancement in the haze region. 

No study that has proposed an astrophysical explana- 
tion for the haze has demonstrated that it can satisfy all 
of these criteria. Pulsars are the primary astro physical solu- 
tion t hat has been discussed i n the literature (IZhang et alj 
120091 ; iKaplinghat et all |2009| ; lHarding fc Abazaiianl l2010h . 
We argue here that a disk population of pulsars likely can- 
not satisfy the third criterion. 

If dark matter annihilations from the Milky Way's hy- 
poth etical dark matter density cusp are the source of the 
haze, iFinkbeinerl (|2004bh found that these three criteria can 
easily be satisfied. Furthermore, for standard dark matter 
models, large clumping factors or boost factors in the an- 
nihilation cross section above theoretically pref erred values 
may not be required to e xplain this radiation (|Finkbeinerl 
l2004bl ; iHooper et ai]|2007l ) . If large boost factors are not re- 
quired, this would be in contrast to other recent astrophys- 
ical anomalies that have been interpreted as arising from 

dark matter annihilations. 

Recently, a 7-ray "haze" was detected by IDobler et al.l 

120101 ) using data from the Fermi Space Telescope, which has 



a similar spatial morphology to the WMAP haze. This haze 
was interpreted as inverse Compton emission off of starlight 
from the same hard popula tion of electrons th at generate 
the haze in the microwave (|Dobler et alj l20icl b The 7-ray 
haze may be another difficulty for astrophysical models of 
the haz e that likely can be accommodated in dark matter 
models l|Cholis et alj2009al l. However, we argue that present 
measurements of this excess in the 7-ray do not yet provide 
a stringent test for models of the haze. 

The bulk of this paper investigates different explana- 
tions for the haze emission and their observational signa- 
tures. Section [5] summarizes the standard diffusive model 
for galactic cosmic ray propagation. Section [3] details sev- 
eral candidate sources for the haze. Section [4] contrasts the 
microwave properties of the candidate haze models, and Sec- 
tion [5] contrasts these models' inverse Compton signatures. 
Throughout this article we refer to the radiating particles as 
"electrons", even though some fraction may be positrons. 



2 GALACTIC COSMIC RAY PROPAGATION 

Cosmic ray electrons travel through the Milky Way halo, 
scattering off of small-scale magnetic inhomogeneities. While 
there is significant uncertainty in the character of cosmic ray 
trajectories, the path of cosmic ray electrons is typically ap- 
proximated as a random walk because the distance a cosmic 
ray electron travels in one direction (approximately the Lar- 
mor radius or 10 13 [E/30 GeV] {B ± /WnG}~ 1 cm) is much 
smaller than the size of the galaxy and its radio halo. In 
this approximation, the electron number density n evolves 
according to 



^ = ^ • (K{E, x) ^n) + A {B (E, x) n) + Q(E, x), (4) 

where K and B are diffusion coefficients, Q is the 
source term, and x r epresents the galactic location (e.g. 
iBerezinskii et al.lll990l )FI 

Under the simplistic assumption that the diffusion pa- 
rameters are spatially independent, the number of cosmic 
ray electrons at distan ce r away from a burst that occurred 
a time At ago is fe.g.. IBerezinskii et al.|[T990l ) 



G(r, I A) 



r \21-3/2 

27rA exp 



r 

~2A2 



(5) 



where the diffusion length A is given by 
[■E a 

A 2 = 2 / dEK(E)B(E)~ 1 , 

J E 

2K [l - (1 - AEAt) 1 -^ 

A(1-S)E 1 ~ S ' 

and E = E/(l — AE At) is the input energy. We have 
adopted the standard parameterizations B(E) = AE 2 (the 



(6) 



1 Equation J2) assumes that the diffusion is isotropic and that 
re-acceleration within the interstellar medium is unimportant. 
Re-acceleration is thought to be u nimportant at E relevant to 
this study (e.g.. I&trong et aT1l2007t) . Cosmic ray data constrains 
re-acceleration of nuclei with an energy per nuclear charge of 
> 1 GeV to be negligible, and re-acceleration should be even 
less rele vant for electrons at the se energies owing to their shorter 
lifetime dBerezinskii et alj[l990h . 
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scaling for inverse Compton and synchrotron losses) and 
K(E) = K (E/lGeV) s , where S > (such that more en- 
ergetic electrons diffuse faster). Also relevant is the electron 
distribution for a burst with an input spectrum of QoE~ liat 
rather than a 5- function at Eo- In this case, the distribution 
of electrons at E < (A At) -1 is 



N(r, E\At) = Q G(r [ A) (1 — AEAt) 



Tint -2 



E~ 



(7) 



In the limit that E ~ Eo (equivalent to AE At <C 1), 
equation @ reduces to 



A = 2.3 



Ko 



10 29 cm 2 s- 1 



/ At 



V 10 6 yr 



— V 

30 GeV J 



kpc. 



(8) 

This limit is most relevant for characterizing the diffusion of 
an event in which all of the energy in cosmic rays was input 
at a single time. In the opposite limit AE At — > 1, equation 
((6]) reduces to 



A = 3.6 



K /A 



10 44 cm 2 GeV 



30 GeV 



kpc. (9) 



This second limit is the relevant limit for a source with a 
continuous output at energy Eo S> E (which approximates 
the input if dark matter annihilations are the source; Section 
13. 5|) . We find the the diffusion length has essentially the 
same scalings as in equation (J5J) for a continuous input of 
cosmic rays with a power-law spectrum. Although, for an 
input power-law of 1.5 < Tint % 3, the coefficient in front of 
equation (O becomes « 1.5 — 2.5 rather than 3.6. 

Both equations ((HJ and are evaluated at 5 = 0.6, and 
the evaluated values for S, Ko, and A are consistent with 
measurements as described below. The diffusion distance in 
both of these equations is comparable to the radial extent 
of the haze or ~ 3 kpc. Thus, it is possible that source(s) in 
the center of the Galaxy could be responsible for the haze. 

A key difference between the continuous source case 
and the case in which a single event was the source is that 
presently less energetic electrons (which were released at an 
earlier time) have diffused further for a continuous source 
(A ~ E (5 ~ l)/2 where 0.3 < S < 0.6). The opposite trend 
with energy holds for an event (A ~ E s ' 2 ). These different 
scalings arise because for an event A is set by how far elec- 
trons traveled in At, whereas for a continuous source A is 
set by the distance electrons of energy E traveled prior to 
losing their last AE ~ E of energy. In addition, the aver- 
age spectral index of the electrons for an event is equal to 
Tint f° r E <C (AAt) -1 , and there are zero electrons with 
E > (AAt) -1 (eqn. [7]). In contrast, for a continuous source 
the average index is equal to Tint + 1 for Tint > 1 and equal 
to 2 otherwise. 

To match cosmic ray observations, galactic cosmic ray 
models assume that cosmic r ays escape freely a t some dis- 
tance Lhalo above the disk (|Strong et al.l [2007). One can 
generalize the solution to the diffusion equation, G(r|A), to 
satisfy such a boundary condition by the method of images 
with a convergent series of images. Furthermore, the cos- 
mic ray distribution for a spatially or temporally extended 
source can be derived via convolution with G(r\X). We adopt 
this approach in most of our calculations. The drawback to 
this approach is that the coefficients of the diffusion must be 
spatially constant. We note that even the most sophisticated 



models for cosmic r ay diffusion, such as the GALPROP code 
jStrong et alj|2007l ). assume that K is independent of x. 

Unless stated otherwise, we use B± = \/2/3 10 /xG 
and A = 10 -15 GeV -1 s -1 . The latter is equivalent to 
U T + Ub ~ lOeVcm -3 and is several times the value at the 
solar circle. This value for B± is consiste nt with the esti- 
mated value at ~ 3 kpc in Galactic radius (Beck 2001). The 
parameter A is degenerate with other diffusion parameters 
(both with K and with At). In addition, we adjust Ko be- 
tween 0.1 and 2 x 10 29 cm 2 s _1 and 8 between 0.3 and 0.6 in 
order to qualitatively reproduce the haze observations in the 
context of different source models. The former value of 0.3 
is what is expected for Kolmogorov spectrum of magneto- 
hydrodynamical turbulence. The latter value of 0.6 is closer 
to what is expected for a Kraichnan-type turbulent spec- 
trum. In addition, we set Lhaio = 4 kpc unless stated other- 
wise. The standard diffusive cosmic ray models are able to 
fit present cosmic ray data with 5 = 0.3 — 0.6, L na io = 3 — 5 
kpc, and Ko = 3 - 5 x 10 28 cm 2 s" 1 ([Strong et al.ll2007l ). 
Although, larger Ko can be accommodated if Lhalo is also 
increased, and Ko could be as mu ch as a factor of 10 larg er 
in the radio halo than in the disk l|Berezinskii et al.l ll990')l 2 l 



3 POTENTIAL HAZE SOURCES 



To detect the haze, first iDobler fc Finkbeinerl (|2008h re- 
moved the Cosmic Microwave Background (CMB) temper- 
ature anisotropics from the WMAP sky maps, which can be 
much larger in amplitude than the foregrounds. They took 
linear combinations of the different frequency channels that 
roughly add up to zero in brightness temperature to perform 
this removal. There is uncertainty in how to do this correctly 
since the optimal linear combination depends on the prop- 
erties of the foregrounds that one is measuring, and this 
uncertainty dominates the measureme nt errors on the haze 
profile . To quantify this uncertainty, IDobler fc Finkbeinerl 
(2008) used 6 physically motivated estimators for the CMB 
(i.e., different linear combinations) to subtract the CMB 
anisotropies. The shaded regions/points for the haze mea- 
surement that are included in some of the ensuing plots 
quantify the range of measured values from these different 
estimators. 

Once the CMB anisotropie s were removed, this yielded 
a map of the foregrounds. Next . IDobler fc Finkbeiner (|200&f ) 
simultaneously fit templates for galactic dust emission, free- 
free (using Ha emission as a tracer) , and "soft" synchrotron 
(using 408 MHz data). The haze was the residual emission 
that these templates did not account for. Thermal dust is not 
a significant contaminant of the lowest WMAP frequency 



2 The constraints on the parameters of this standard diffusive 
model primarily derive from two measurements. First, measure- 
ments of the ratio of different primary and secondary cosmic rays 
yield the column of protons that cosmic rays traverse (which re- 
lates to the distance traveled). Second, the abundance of radioac- 
tive cosmic ray nuclei provide the lifetime of cosmic rays in the 
galactic diffusion zone. However, it is likely that the cosmic ray 
diffusion parameters are different above the disk in the central 
< 4 kpc of the Galaxy from their effective value obtained from 
local measurements, but even the sign of how Ko (for example) 
will differ from its locally-measured value is uncertain. 
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channels at which the haze detectio n is most secure (al- 
thoug h spinning dust can contribute; iDobler fc Finkbeinerl 
2008), and Ha is a robust tracer of free- free emission from 
~ 10 4 K gas (provided that absorption by intervening dust 
is minimal). In addition, direct and indirect emission from 
unresolved sources in the g alactic bulge is also unlikely 
( Ban dvopadhvav et al.l [20 08). Thus, the most compelling 
models for the haze are either that it is free-free from gas 
with T S> 10 4 K or that it is additional synchrotron radi- 
ation that is not c aptured in the soft-synchrotron template 
l|Finkbeinerll2004al ). This section discusses candidate sources 
for the haze that fall into these two categories. 



3.1 Free-Free Emission 

Soft X-ray maps of the inner g alaxy show a similar spatial 
morphology to that of the haze l|Snowden et al.|[l997l ). which 
is suggestive that the haze owes to free-free emission from 
hot gas. However, previous studies have shown t hat it is dif- 
ficult to create the haze with free- free emission (|Finkbeinerl 
l2004bl : IDobler fc Finkbeinerl [2008). In fact, if the haze a lso 
appears in the 7-ray, as suggested bv lDobler et alJ ([2010T ). it 
cannot be free-free. We review the arguments for why this 
microwave emission most likely cannot arise from free-free 
processes. 

Firstly, IDobler fc Finkbeinerl (2008) claimed that the 
radio spectral index of the haze (defined here as the power- 
law slope of the specific intensity) is inconsistent with free- 
free emission. The radio spectral index o f free-f ree radiation 
is ~ 0.15, whereas IDobler fc Finkbeinerl (|2008l ) measured a 
median radio spectral index for the haze of ~ 0.4 between 
the 23 or 33 and 4 1 GHz channels. Howe v er, th ere were still 
estimators in the IDobler fc Finkbeinerl (|2008l ) study that 
yielded a spectrum that was hard enough to be consistent 
with free-free emission. Even if their claim is correct that the 
emission is too soft to be free-free, a component of the haze 
could arise from free-free emission, and this component plus 
a softer component could yield the observed haze spectrum. 

If a significant component of the haze were due to free- 
free emission from < 10 5 K gas, the sam e gas would be vis- 
ible in Ha radiation l|Finkbeinerll2004al ). and, if it were due 
to > 10 6 K gas, it would be visi bl e in the ROSAT so ft X-ray 
maps of ISnowden et all (|l997l ) (|Finkbeinerll2004bl ). Thus, 
only free-free emission from ~ 10 5 -10 6 K gas is viable. How- 
ever , gas at suc h temperatures is typically thermally unsta- 
ble (Field 1965), and to produce the haze emission requires 
a large column density of ions at these temperatures (1 kpc 
with a density of 0.1cm -3 ). 

Energetically it is difficult for there to be a sufficient 
amount of gas in this range of temperatures to produce the 
haze. Since the free-free specific intensity is approximately 
independent of frequency for hu < kt,T g (where T g is the 
gas temperature), the free- free energy loss rate of such gas 
would be 

E blBm ~ K„] taj ^- ~ 10 56 (-^-) ergMyr- 1 , 

h fhaze \ 10 b K J 

(10) 

which requires the kinetic power of ~ 10 5 supernovae ex- 
plosions per Myr for T g = 10 6 K gas (and larger powers 
are required if resonance line emission is also included). It 



is unlikely that Galactic Center sources can supply these 
energetics. 

Figure Q] illustrates the constraints in the gas temper- 
ature (T g ) versus emission measure (EM) plane if free- free 
emission were the source of the haze, where EM is defined as 
the line-of-sight integral of the square of the electron num- 
ber density. The haze must fall on the solid line in this plane 
to yield the measured amplitude at b = 10° . The shaded re- 
gion at the left is excluded by Ha measurements and at the 
upper-right by soft X-ray measurements^ The long dashed 
curves are contours of constant energy loss for gas at the 
solar metallicity, with the curve just above the solid line 
representing 10 56 erg Myr -1 , and the other curves step in 
factors of 10 from this value. In regions rightward of the dot- 
dashed vertical line, the gas cannot be gravitationally con- 
fined in the Galactic Plane, and, therefore, this parameter 
space is also excluded. See the figure caption for more de- 
tails. Because of the required energetics, gas temperatures, 
and pressures, it is unlikely that the haze emission is due to 
free- free. 



3.2 Synchrotron Radiation from Cosmic Rays 
Produced by the Disk Population of 
Supernovae 

IDobler fc Finkbeinerl <|200Sl ) used the Haslam 408 MHz sky 
map as a template to subtract the "soft" synchrotron com- 
ponent of the foreground emission. This procedure assumed 
that the distribution of ~ 3 GeV electrons that is responsi- 
ble for the 408 MHz emission is the same as the distribution 
of ~ 20 GeV that is responsible for the soft synchrotron 
emission in the WMAP bands. Even if the sources of the 
electrons at both energies are the same, this assumption 
can be violated by diffusion effects: The diffusion length for 
a 20 GeV electron is 1.4 — 1.9 times smaller than a 3 GeV 
electron for 0.3 < 5 < 0.6. 

The IDobler fc Finkbeinerl (|200Sl ) analysis found that 
the soft synchrotron emission at 23 GHz is roughly 
{1,2,3} kJy sr -1 at {25,10,5} de grees in latitude and 
at longitude / = 0, whereas the haze is measured to be 
{1,2,5} kJy sr^ 1 . At fixed angle from the Galactic Center, 
but for I 7^ 0, the soft synchrotron fraction is larger. To cre- 
ate the correct amplitude for the haze at / = 0, a hardening 
of ~ 0.25 in the synchrotron spectral index (~ 0.5 in 7) is 
required when extrapolating from 408 MHz. 

It would be difficult for diffusion effects from the disk 
population of supernovae to result in a change in 7 of 0.5 
in the inner 20° . The distribution of galactic supernovae is 
measured to be a much smoother function of galactic lon- 
gitude than the haze emission, with the favored model pre- 
ferring a distribution that peaks at 4 kpc or I = 30°, and 
wit h few supernovae wi thin 2 kpc from the Galactic Cen- 
ter jLorimer et al.l 120061 1 . The supernova distribution is de- 
rived from emission measure determinations from hundreds 
of galactic pulsars and a model for the Milky Way interstel- 
lar electron density. 



3 This region only accounts for the contribution of free-free emis- 
sion to the X-ray flux. The inclusion of the not insignificant emis- 
sion from metal lines would shift this contour down and to the 
left, reducing further the allowed region. 
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Figure 1. Red solid curve represents the required T g and EM if 
the haze owes to free- free emission (assuming I v = 2 kjy sr~ , 
the haze value at b = 10°). The gray region on the left is excluded 
by Ha measurements (with the limit > 0.3 Ralcighs; Fin kbeinerl 
l2004al) and on the right by X-ra y measurements (> 7 Jy sr _1 
at 750 eV; ISnowden et al.l Il997t) . The dashed curves are con- 
tours of constant energy input E, using the approximate formula 
E = 2ir /3 A EM A;j,„„, where A is the co oling rate density at solar 
metallicity (Sutherland & Dopita 1993) and we set Aj lazc = 2 kpc 
— the approximate extent of the haze. The dashed curve that 
falls just above the solid red curve is for 10 56 erg Myr -1 , and the 
other dashed curves represent E that are multiples of 10 of this 
value. The dot-dashed vertical line represents the temperature 
that balances the radial gravity of the disk for a circular velocity 
of v c = 100 km s~ x (roughly the Milky Way circular velocity at 
2 kpc) derived using kT g = 1.2 m p v^. /2. 



However, the inferred distribution of pulsars depends 
strongly on modeling of the interstellar electron distribution. 
lLorimer et al.l (|2006l ) showed that a more centrally concen- 
trated distribution of pulsars can be obtained by assuming a 
less physically motivated, strongly peaked electron distribu- 
tion towards the Galacti c Center. Yet, even in the most ex- 
treme case considered in lLorimer et all (|2006l ). the inferred 
radial distribution is relatively flat in the central < 5 kpc, 
with ~ 3% of disk pulsars in the inner 1 kpc. 

For a source distribution that peaks 4 kpc into the disk, 
or even a flat radial distribution in the inner several kpc, 
most of the galactic cosmic ray electrons must diffuse from 
their production sights a few kpc to reach the Galactic Cen- 
ter region. Since less energetic electrons will have diffused 
further on average, diffusion causes the synchrotron emission 
from this population to soften toward the Galactic Center. 
We used the GALPROP code to test this argument, as well 
as to test how well the 408 MHz map can be used as template 
for the microwave synchrotron from supernova-accelerated 
electrons. In part icular, we used the d efault GALPROP pa- 
rameters given in lStrong et al.l (|2004l ). as well as slight vari- 
ants of these, to generate mock synchrotron sky maps. GAL- 
PROP is a publicly available code that provides a numerical 
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Figure 2. Residuals after subtraction of the fiducial GALPROP 
model prediction at 408 MHz from its prediction at 23 GHz, where 
the 408 MHz prediction has been renormalized to minimize the 
square of the map differences. The solid curve is for I = 0, the 
dashed for I = 25°, and the dotted for I = 40°. This procedure 
results in a small negative residual towards the Galactic Center 
at latitudes b ~ 10° for the standard distribution of cosmic ray 
sources. Its amplitude is significantly smaller than the amplitude 
of the haze. 

solution to equation |4) in a cylindrical region The de- 
fault parameters are tuned to fit many different cosmic ray 
observations, using for the sources the standardly assumed 
distribution of supernovae. In particular, the default model 
takes a radial distribution of supernovae in the disk that 
scales as R 4 exp(— R/S kpc). 

We subtracted the 408 MHz map generated with GAL- 
PROP from the 23 GHz GALPROP map, minimizing the 
square of their difference in the region defined by —45 < 
I < 45° and < b < 4 5° to emulate the method in 
iDobler fc Finkbeinerl (|2008l ). The results are shown in Fig- 
ure [5] This procedure resulted in a small negative signal 
toward the Galactic Center region at latitude b ~ 10° and 
longitude I — (solid curve, Fig. [2}. This signal becomes less 
negative with increasing longitude (the dashed and dotted 
curves show I = 25° and / = 40°). This residual is an order 
of magnitude smaller in amplitude than the haze0 In addi- 
tion, we find an almost identical result if we instead assume 
a flat radial distribution in the inner several kpc of the disk, 
with source density oc exp (— R/8 kpc), a n alogo us to the ex- 
treme case considered in lLorimer et al.l i|2006l ). While the 
GALPROP cosmic ray propagation model is simplistic, this 
small negative signal stems from the distribution of sources 
that is assumed and should be robust to the propagation 
model. 

In general, a disk-like population of sources does not 
produce a spectrum where the spectral index varies signifi- 
cantly across the sky as long as the magnetic field is a weak 

4 http:/ /galprop. stanford.edu/web_galprop/galprop_home. html 

5 The divergent component at small angles in Figure[2]owes to the 
fact that, for a continuous source population, the number density 
of electrons diverges as the distance to the source (which here 
is the disk of supernovae) decreases. This divergence in Figure [2] 
affects smaller angles than are relevant for the haze emission and 
should be ignored for this discussion. 
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function of radius and the solar circle is far from the radial 
edge of the diffusion region. (These criteria appear to be 
roughly satisfied in the bulk of the Galaxy owing to the ap- 
proximately constant spectral index that is observed across 
the bulk of the sky between the 408 MHz maps and those in 
the microwave.) To understand this, imagine different cylin- 
drically symmetric diffusion zones that characterize different 
regimes for the electrons/magnetic fields around the disk of 
sources. All sightlines from within this disk (that do not 
intersect the radial boundary) will penetrate each of these 
diffusion zones in the same proportion, resulting in the same 
spectral shape for the electron (and synchrotron) spectrum 
as a function of angle. 

In contrast, for a point cosmic ray source in the Galactic 
Center, the diffusion zones are more spherical, and sightlines 
will intersect these zones in different proportions. The rest 
of this section investigates different point source models for 
the cosmic ray injection. 

3.3 Synchrotron Radiation from Cosmic Rays 
Produced by Supernovae in the Galactic 
Center 

We argued that the supernova population inferred from ra- 
dio pulsar surveys likely cannot be the source of the haze 
emission. However, there is a population of undiscovered 
pulsars in the Galactic Center (indicating additional super- 
nova activity) that has been missed by radio pulsar surveys 
and that are not incl uded in galactic cosmic ray models. 
ICordes fc Laziol l|l997h demonstrated that previous pulsar 
surveys would not have identified pulsars in the central cou- 
ple hundred parsecs because pulse broadening from dense 
gas in the Galactic Center washes out the pulses at the sur- 
veyed frequencies. 

In addition, the observed star formation rate in the 
Galactic Center can produce enough supernovae to satisfy 
the haze energetics. For Ub + U r = 10 eV cm -3 , 1000 su- 
pernovae in the Galactic Center must have occurred in the 
last 10 6 yr in order to source the required ~ 10 51 erg in elec- 
trons if ~ 0.1% of the kinetic energy of each supernova is 
converted into GeV electrons that reach the ISM. This per- 
centile of 10 51 erg (or 10 48 erg) is the amount of energy that 
has been estimated to go into E > 1 GeV electrons from 
SN 1006 and also from observations of the intergalactic cos- 
mic ray electron spectrum comb ined with estimates for the 
Milky Way supernova rate (see iKobavashi et al . 2004 and 
references therein) One thousand is ~ 10% of the number 
of supernovae that occurred in the Galaxy within this pe- 
riod, which is a reasonable fraction since 10% of the present 
galactic star formation is esti mated to re side in the inner 
couple hundred parsec region (Figcr 2008). 

A complementary estimate for the number of required 
supernovae can be derived by noting that the haze is com- 
parable in amplitude to the emission from soft synchrotron 
(which is attributed to supernovae) within the inner 20° , and 
soft synchrotron emission of comparable amplitude covers 
/sky of the sky, where / a k y ~ 0.1. Galactic electrons at these 

6 I Thompson et aL I j2006f ) derived a largely model-independent 
estimate of f» 10 49 erg per supernova from the far infrared-radio 
correlation in starburst galaxies. 
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Figure 3. Supernova models for the WMAP haze. Top panel: 
The intensity profile as a function of angle from the Galactic 
Center for two supernova models described in the text (dot- 
dashed and dashed curves) as well as the measured intensity pro- 
file of the WMAP haze in the 23 GHz channel (solid black curve 
with dots, with the grey region quantifying the current measure- 
ment uncertainty). The thin curves are our models' raw intensity 
at 23 GHz. The thick curves are this minus their intensity at 
408 MHz scaled to 23 GHz, assuming I u ^ v ^ to rescale the 
408 MHz map, emulating t he template subtraction procedure in 
iDobler &: Finkbeinen 1 120081) . Bottom panel: The ratio of the in- 
tensity in the 23 and 33 GHz channels. The curves show the same 
models as in the top panel. The dotted horizontal curve is this 
ratio for l v ~ v~ 1 — indicative of the galactic soft-synchrotron. 
The filled circles represent the hardness values as quantified by the 
ratio of the total haze intensity in the 23 and 33 GHz channels for 
the 6 different estimators to remove the primordial anisotropics in 
IDobler &: Finkbeinerl J2008h ■ The leftmost column of circles uses 
the power-law index inferred between the 23 and 33 GHz chan- 
nels, the middle that between the 33 and 41, and the rightmost 
that between the 23 and 41. 

energies are expected to lose roughly half of their energy by 
synchrotron emission (and roughly half by inverse Comp- 
ton). Therefore, the fraction of Galactic supernovae required 
to produce the haze is ~ (20°) 2 /[4ttx (180/vr) 2 x/ sky ] w 0.10 
If supernovae are the source of the haze, diffusion ef- 
fects are required to harden the electron spectrum in the 
haze region. Galactic Center supernovae can be treated as 
a continuous source population such that electrons that 
have diffused the furthest will have cooled the most (eqn. 
[9]). Therefore, diffusion results in the synchrotron spectrum 

7 These estimates assume that ~ 30 GeV electrons can escape 
from the central couple hundred parsec regi on. There are som e 
indications of milli-Gauss fields in this region ( Plant e et al Il995l h 
and, if this were the case, haze electrons would likely not be able 
to escape from this region. However, it is probable that milli- 
Gauss fields that have been detected are localized to bundles and 
that in most of the volum e the magnetic fields are comparable to 
in the rest of the Galaxy llUchida fe Guesten|[l995h . 
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from Galactic Center supernovae softening as the distance 
from the Galactic Center increases. At r < A, the haze elec- 
trons should be hardened in this model by 0.6 < A7 < 
1.1 for 0.3 < S < 0.6 noting that A ~ E {s ~ 1)/2 and 
n(E', r) /n(E, r) ~ \(E) 3 / X(E') 3 as small r. The thin curves 
in the left panel in Figure [3] show the intensity profile for 
two simple supernova models; dashed: Kq — 2 x 10 29 cm 2 
s _1 , 7i nt — 2, S — 0.3, and with input electron cosmic 
ray luminosity E to t = 2 x 10 38 erg s _1 (equivalent to 
6 x 10 51 erg Myr" 1 ); dot-dashed: K = 1 x 10 29 cm 2 s~\ 
7int = 2,6 = 0.6, and E tot = 3 x 10 38 erg s" 1 . To calculate 
i5tot, the input spectrum of cosmic rays is assumed to extend 
between 1 GeV and 1 TeV. Both models take a value for Ko 
that is slightly higher than is typically used in cosmic ray 
diffusion calculations. (Although, A might also be on the 
high side, and the quantity of importance for diffusion is 
Kq/A.) The choice of 7mt = 2 for supernovae is empirically 
motivated^ 

Both of these models are normalized to over-predict 
the measured haze intensity (the solid black curve with 
dots in the top panel in Fig. [3| for reasons detailed be- 
low. Both models also predict a softer spectrum as quan- 
tified by the ratio of the WMAP intensity in the 23 and 
33 GHz channels than the inferred estimates for this ra- 
tio from lDobler fc Finkbeineil (|2008f ). The inferred estimates 
use the measured haze 7 between three permutations of the 
lowest three WMAP frequency channels plus six different es- 
timators to remove the CMB anisotropies (the filled circles 
in the bottom panel of Fig. [3j> 13 

The thick curves in Figure[3]may be closer to what is ac- 
tually measured than the thin curves. These curves represent 
the intensity after subtracting the specific intensity I v of the 
corresponding model at 408 MHz scaled to 23 GHz assuming 
In oc v ^ , the measured average radio spectral index for the 
Galactic soft synchrotron. This subtraction should be simi- 
lar to how the haze is actually obtained. The thick curves in 

8 Numerous observations have been made of the galactic popu- 
lation of electron cosmic rays at GeV energies, which are believed 
to be accelera ted in supernova re mnants. The fiducial cosmic ray 
model used in lStrong"e t al. (2004) and in most GALPROP calcu- 
lations assumes that the injected supernova spectrum transitions 
over ^2 — 20 GeV from 7; nt = 1.5 at lower energies to 7; nt = 2.5 
at higher energies, with an effective injection spectrum between 
3 and 20 GeV of 7i nt ~ 2. This extended transition is used to 
fit a break in the lo cally measured cosmic ray electron spectrum 
llStrong et al,l [20041) . (Although, at least some of this observed 
break may be explained by the increased importance of collisional 
losses at a couple GeV.) Some cosmic ray models have assumed 
a single power-law of 7i nt = 1.9 across all re levant energies, an d 
these models are able to fit much of the data l lStrong et ay 2000). 
Empirically, the spectral index of the soft synchrotron between 
the radio and microwave (and across different microwave bands) 
is measured to be 7 Ri 3 on average, the anticipated spectrum 
after cooling if 7; nt = 2. 

9 Unfortunately, present measurements of the haze are not able 
to constrain the hardness as a function of radius, which would 
be the ideal discriminator of the models presented in this pa- 
per. What has been measured is the amplitude of the emission in 
the different frequency channels assuming a template with profile 
I u {6) ~ 1/0 - 1/0 O for 5 < < O = 45°. The average intensity- 
weighted angle for this profile is (0) = 18°, which is the angle at 
which we plot the intensity ratios in the bottom panel of Figure 
® 



the top panel in Figure [3] are able to largely reproduce the 
profile of the haze emission at 23 GHz for 9 < 20° . However, 
these haze models produce a spectral hardness that is only 
consistent with the couple estimators that yield the softest 
values for the haze hardness (bottom panel, Fig. |3J. 

There are two additional potential problems with the 
Galactic Center supernova explanation for the haze. One 
problem is that to reproduce the amplitude of the haze at 
23 GHz requires an amplitude prior to subtracting the soft 
synchrotron that is 150 — 200% of the amplitude of the haze 
in our two models in Figure [3] (thick curves). A value of 
200% is close to the total observed emission in the haze 
region at 23 GHz, and such emission will significantly cut 
into the budget for the soft-synchrotron contribution from 
the disk population of supernovae. The second problem is 
that the electrons in our two supernova models contribute 
to the emission at 408 MHz at the 25 — 50% level in the haze 
region. However, there is no significant centralized emission 
in the inner ~ 20° in the 408 MHz maps, at least at the 50% 
level. 

One way to alleviate both of these potential issues is if 
the diffusion length were longer than we have assumed such 
that the haze spectrum from supernovae were even more 
hardened by diffusion (Section [3J. Another way is if there 
were a break in the injected electron spectrum of supernova 
between 5 GHz and 20 GHz such that the spectrum is harder 
at lower energies. Such a break may even be required to 
fit local cosmic ray data. 9 Such a break could increase the 
average diffusion length of electrons at 408 MHz, making 
the emission in this band from the Galactic Center sources 
more diffuse. However, this solution would also reduce the 
additional hardening our supernova models gain during the 
subtraction of the soft synchrotron. 

3.3.1 Pulsars Instead? 

Most previous astrop hysical explanation s for the haze have 
focus ed on pulsars JZhang et all 120091 : iKaplinghat et al.l 
2009: lHarding fc Abazaiianll20ld ). While it is not clear how 
efficiently pulsars can produce interstellar cosmic ray elec- 
trons, several studies have attributed recently discovered 
cosmic ray "anomalie s" in the ATIC JChang et al.l 120081 ) 
and Pamela data set s jAdriani et a l. 2009) to these objects 
(e.g., IProfumol [20081 : iHooper et al.ll2009l ). This led to sev- 
eral studies investigating whether pulsars could also explain 
the WMAP haze anomaly. However, the generation of the 
haze with pulsars may be incompatible with the explana- 
tion that the ATI C and Pamela anomalies owe to pulsars 
([Zhang et al.ll2009l ). 

Young pulsars (< 10 4 yr) have been conceived as a 
promising model for the haze because they are conjectured 
to supply a harder source of electrons and positron s at 
these energies than supernov ae. IKaplinghat et al.l (|2009l ) and 
lHarding fc Abazaiianl (|2010j ) argued that a disk population 
of young pulsars could be responsible for the haze if these 
pulsars output ~ 10% of their spin-down energy into inter- 
stellar cosmic ray electro ns and have 7i n t ~ 1-5. In addition, 
IKaplinghat et all {2009) showed that varying the longitu- 
dinal distribution of the uncertain galactic magnetic field 
could produce a similar intensity profile to the haze even 
from a disk population of pulsars. 

However, previous analyzes that invoked pulsars to ex- 
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plain the haze largely neglected that, when the soft syn- 
chrotron is subtracted from the WMAP maps using the 
408 MHz data as a template, this procedure should remove 
much of the pulsar contribution in the microwave because 
young pulsars have essentially the same spatial distribution 
as supernovae F°l Therefore, pulsar models also require a dif- 
fusion hardened electron spectrum in the Galactic Center to 
produce a haze-like synchrotron signature. The same argu- 
ments in Section [372] for why a disk population is unlikely to 
produce the haze also apply in this case. 

3.3.2 Supernovae with a Galactic Wind 

A wind would naturally harden the spectrum towards the 
Galactic Center by transporting cosmic rays that have lost 
the most energy the furthest In the limit of zero diffusion 
and zero adiabatic losses, the electron spectrum of a wind 
should have a spectral index of 7 = 7i nt for E < (Ar/v)^ 1 
and have a cutoff due to cooling above this energy, where r 
is the distance from the Galactic Center and v is the wind 
velocity. For Ub + U r — 10 eV cm -3 , a wind can harden 
the electron spectrum over scales comparable to that of the 
haze if it can transport particles > 2 kpc within 1 Myr (the 
timescales for the electrons that would produce the haze to 
lose their energy) . To meeting these conditions would require 
a 2000 km s" 1 wind. 

There may be evidence for a galactic wind in the Milk y 
Way (iBland-Hawthorn fc Cohen|[20ol ; lEverett et aT1l2008h . 
lEverett et alj (^OOct ) interpreted a ROSAT soft-X-ray excess 
with similar extent to the haze as indicating a wind, and 
their best fit models yield v — 200 — 300 km s _1 . However, 
for v = 300 km s" 1 and U B + U r = 10 eV cm" 3 , 30 GeV 
electrons could only reach a radius of 0.3 kpc — much less 
than the extent of the haze. A haze-like extent would require 
a significantly smaller value for Ub + U t and/or a much faster 
v. In addition, these estimates ignore adiabatic losses, which 
would be present for any expanding outflow. We conclude 
that it is unlikely that a wind alone can produce the haze 
(unless there is in-situ electron acceleration in the wind). 

3.3.3 Other Source Populations in the Galactic Center 

The discussion in subsection 3.3 also applies to other pop- 
ulations of sources in the central ~ 1 kpc of the Galaxy. 
Any such source population would have the same difficulty 
as supernovae of being consistent with measurements of the 

10 iKap linghat ct al. (2009) suggested one mechanism that could 
decouple the spatial distribution in these bands: If the domi- 
nant contribution to the emission at 408 MHz owes to cosmic 
ray secondaries rather than primaries. Local measurements of the 
positron fraction suggest that ~ 10% of the e+e- that contribute 
at 408 MHz owe to secondaries. Therefore, while possible since 
local measurements may not be representative of the galaxy, this 
hypothesis is not supported by present data. It is also unclear 
how different the distribution of secondary sources, which should 
largely trace the gas distribution, will be from that of the pri- 
maries, which traces star formation (which to some extent traces 
the gas). 

11 The haze electrons themselves cannot drive a strong wind. The 
pressure from the haze electrons is minute compared to that from 
the much more abundant cosmic ray protons. 



haze hardness unless they had an injection spectrum with 
Tint ;$ 2. For example, tidal disruptions by stars that fall 
close enough to the supermassive black hole are estimated 
to occur every 10~ 4 to 1 0~ 5 yr~\ with kinetic en ergies 
of to 51 " 52 erg per event jStrubbe fc QuataertJ [jpoj ). If a 
small fraction (~ 1%) of this energy goes into accelerating 
electrons, the energetics would be sufficient to produce the 
haze. Secondly, the accretion rate onto Sagittarius A* is esti- 
mated to be a few x 1 ~ 6 from wind s off of nearby stars 
|Quataert et al.lll999l ; iGenzel et ai1ll994h . and upper limits 
on this rate based on observations coupled with accretion 
mode ls are 1 — 2 orders of magnitude larger (|Quataert et al.l 
Il999l) . If even a small fraction of this rest mass energy con- 
tributed to accelerating electrons (possibly via an associated 
jet), this low level of accretion would also be sufficient to 
power the haze. 

3.4 A Single Energetic Event in the Galactic 
Center 

The next possibility that we consider is whether a single 
explosive event could be responsible for the haze. Such a 
transient must output > 10 5 erg (eqn. 1; and we find ~ 
10 52 erg in Section |3} . For this model to be different from 
that of a continuous source population, this event must be 
sufficiently rare in order for a cooling break — a chromatic 
break in the electron spectrum to a power-law index of 7i nt + 
1 from 7i nt — to not appear in the microwave from the 
accumulation of electrons from previous events. The absence 
of a cooling break at relevant energies also enables such an 
event to have a hard enough spectrum even with 7i nt > 2, 
in contrast to continuous models. 

For electrons to be present that radiate significant en- 
ergy at frequency v ma m and for a cooling break to not have 
developed at v > v m i n , such a blast must have occurred a 
time 

. 1 _ f B± 50GHz\ 1/2 ( Ur + U B V 1 
Tbiast < 1-0 — — — — 33 Myr 

V 10 flG fmax / V 10 eV C1T1 / 

(11) 

ago. In addition, the last prior energetic event happened 
t > 1.5 ( ^m^™ ) ^ ma ^ min y /2 Myr (12) 

ago, were Tbi as t,max is the value that saturates inequality 
<|11|) . Events that occur every ~ 1 Myr will often accommo- 
date both inequalities. 

Electrons from an event that occurred 1 Myr ago can 
diffuse several kpc (eqn. [8]), on the order of the radial extent 
of the haze. The most energetic electrons diffuse the furthest 
for an event since A ~ E 6/l2 , such that the electron spectrum 
will harden with radius (eqn. [8]). In addition, such an event 
will likely not be visible in the Haslam maps because there 
are far fewer electrons at a few GeV if the explosion has 
7int < 3 than the cumulative number of electrons from su- 
pernovae (for which 7 ~ 3). 

Possibilities to create the haze from a single event in- 
clude episodes that mark the deaths of massive stars (such 
as a rare supernova) and activity by the supermassive black 
hole (SMBH) in the Galactic Center. However, such an event 
would have to produce ~ 10 52 erg in ~ 30 GeV cosmic ray 
electrons (Section [4]), which is likely too large to be from a 



© 2009 RAS, MNRAS 000, 1 



Testing the Dark Matter Annihilation Model for the WMAP Haze 9 



single stellar event. Whereas, transient activity associated 
with the SMBH in the Galactic Center can more easily ac- 
commodate these energetics. In fact, there is evidence for 
recent activity in the Galactic Center. In the 10 4 Mq Arches 
cluster, which sits 30 pc from Sagittarius A*, massive stars 
that have ages < 2 Myr are noticeably absent, indicating 
that this cluster formed w 2 Myr ago. In addition, two other 
young clusters near Sagittarius A* are estimated to have 
twice this age, one of t hese being the Central Cluster in the 
inner 1 pc (|Figerll200Sl ). It is conceivable that the formation 
of one of these clusters was accompanied by a short episode 
of accretion onto Sagittarius A* . If the 4 x 10 6 M SMBH 
shined at less than the Eddington luminosity during a period 
of gas accretion, this period must have lasted > 0.6/ -1 yr. 
The parameter / is the ratio of the energy that goes into 
accelerating cosmic ray electrons, taken to be 10 52 erg, to 
the energy in electromagnetic radiation. 

One compelling candidate for sourcing such an occur- 
rence is the tidal disruption of a small molecular cloud by 
Sagittarius- A* . Molecular clouds are estimated to be tidally 
disru pted by Sagitt arius A* every ~ 10 7 yr in the Milky 
Way (Sanders 1998). In fact, an explanation for the origin 
of the S-stars around Sagittarius-A* is that they wer e formed 
sever al Myr ago from the debris of such an event l|Sanderd 
1998). The timescale for the activity associated with the 
tidal disruption of a molecular cloud by Sagittarius-A* is the 
time for the cloud to cross the black hole. For a cloud size 
of 0.1 pc and a relative velocity of 10 km s _1 with respect 
to Sagittarius A*, the cloud crosses the SMBH in 10 Myr, 
which is much less than 1 Myr as required. Secular activity 
of this sort may be responsible for the nuc lear activity in 
Seyfert galaxies |Hopkins fc Hernquisdl2006r ). However, it is 
also possible that much of the accretion associated with the 
tidal disruption of a small molecular would occur at small 
enough accretion rates such that the accretion disk would 
be in a radiatively inefficient state. In this case, such phe- 
nomena would be more difficult to detect in other galaxies. 



3.5 Dark Matter 

Cosmic rays generated as the byproducts of dark matter an- 
nihilations is the most exciting proposal for the source of the 
haze. With well-motivated assumptions for the distribution 
of dark matter and its properties, the total luminosity in 
dark matter annihilations from the Milky Way halo can be 
estimated: 



L ann (E) » 2 x 10 37 B (^) 



2 -° Mmw 100 GeV _j 
1O 12 M Mdm &8 S ' 



(13) 

Thi s expression assume s an NFW dark matter halo pro- 
file IjNavarro et al.l l 19961 with total dark matter mass Mmw 
and concentration parameter The scaling for c in equa- 
tion (|13[) is approximate and most valid around c = 12, the 
value that matche s rotation curve data for Mmw = 10 12 Mq 
l|Besla et al.ll2007h . M M w = 10 12 M is on the low side of 
Milky Way mass estimates, with most est imates resulting 
in a value that is a factor of 2 larger (e.g. ISakamoto et al.l 



12 We assume a flat ACDM cosmology with Sljf = 0.27 and 
h = 0.7 to calculate the 2 = virial radius in spherical collapse, 
which is necessary to define c. 



120031 ; iLi fc White! I2OO8T ). M DM is the dark matter particle 
mass, and popular theories for weakly interacting massive 
particle (WIMP) dark matter predict M D m ~ 0.1 - 1 TeV. 
The parameter B is a boost (or suppression) factor that 
may arise from additional dumpiness on top of the smooth 
NFW profile or if the value of the WIMP velocity-averaged 
annihilation cross section (a v) is different in the Milky Way 
than at freeze out (where (crv) ~ 3 x 10 -26 cm 3 s _1 is 
required to match the relic abundance). Recent numeri- 
cal work argues that the contribution to B from clump- 
ing is near unity inward of the solar circle (|Springel et al.l 
l2008l : lKamionkowski et alJ2010h . and standard WIMP mod- 
els typically predict that the Galactic value of {av) will be 
comparable to its value at freeze out. However, numerical 
simulations are far from resolving the smallest dark matter 
structures, and dark matter models with B> 1 have been 
invoked by several recent stud ies to explain several percei ved 
astrophysical anomalies (e.g.. lArkani-Hamed et aDl2009l ). It 
is interesting that for standard values of the halo and dark 
matter properties, the dark matter annihilation luminosity 
(eqn. I13p is comparable to the ~ 10 37 erg s _1 required to 
produce the haze. 

The synchrotron haze emission from dark matter anni- 
hilations depends on the parameterization of the Galaxy's 
dark matter profile, magnetic field structure, and distribu- 
tion of interstellar radiation. Because of these uncertainties, 
adding complications in the dark matter sector of our model 
is not warranted, and our dark matter models simply as- 
sume direct annihilations into pairs (a 5-function injection 
at E = Mdm/2). Such models results in the average spec- 
trum of all electrons having 7 = 2 for E < Mdm/2. More 
physically motivated WIMP models (in which e+e- are pro- 
duced after a particle cascade) generally yield electron spec- 
tra that are slightly softer, and only a fr action (~ 1/3) of the 
annihilation energy goes into leptons (|Hooper et al.l 120071 : 
ICumberbatch et alj |2009). All dark matter models produce 
7 = 2 on average for 1 GeV < E <C Mdm- 

For our fiducial Milky Way halo model, the density pro- 
file scales as r _1 out to « 10 kpc, before steepening. As a 
result (assuming dB/dr — 0), an equal amount of annihila- 
tion power per unit radius is emitted at all inner radii, and 
half of the annihilation power originates from within 5 kpc 
for the fiducial halo model. The top panel in Figure [3] shows 
the cumulative energy output from annihilations as a func- 
tion of radius for four plausible halo profiles. The bottom 
panel in Figure [3] shows what this profile would look like if 
the diffusion length were much less than the radial size of 
the annihilation profile and for a uniform B±. (This profile 
also corresponds to the intensity profile for direct annihi- 
lations into 7-rays.) This figure suggests that the angular 
properties of the annihilation profile depend weakly on the 
dark matter halo model for plausible models (assuming that 
substructure is unimportant). 

Figure [5] shows the haze radial intensity profile for four 
illustrative dark matter models in the 23 GHz channel (left 
panel) and also the haze hardness (right panel), expressed 
as the ratio of the intensity in the 23 GHz and 33 GHz 
channels. These models have been tuned to fit the observed 
radial profile, all with {av) = 3 x 10~ 26 cm 3 s" 1 . The red 
solid curve assumes Mdm = 100 GeV, Ko = 10 28 cm 2 s _1 , 
S = 0.3, and B — 10. The green dot-dashed curve is the same 
as the solid curve but with an inner power-law slope of —1.2 
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Figure 4. Top Panel: Cumulative dark matter annihilation lu- 
minosity in units of the total (assuming our fiducial parameters 
of M M w = 10 12 M , c = 12, and an NFW halo profile, un- 
less stated otherwise). The top curve labelled "adiabatic" instead 
assumes an inner power-law of —1.2. The second curve down in- 
stead assumes c = 20. The third down is our fiducial model, and 
the bottom curve instead assumes Mum = 2 x 10 12 Mq. Bottom 
panel: The angular profile of the intensity in annihilations for the 
same four models. 



(what simulations typically find that include gas cooling; 
iTissera et alEoIfr . The blue short-dashed curve is the same 
as the solid curve but for a 10 times more massive WIMP 
and B = 800, and the magenta long-dashed curve is the 
same as the solid curve but for a value for Ko that is 3 times 
larger. All models produce similar intensity profiles, and the 
hardness of the radiation does not depend strongly on angle 
because significant annihilations occur even at several kpc 
from the Galactic Center. In detail, the less massive WIMP 
models produce a slightly softer spectrum, and all WIMP 
models are hardest within r < A from the Galactic Center. 

It is simple to understand more quantitatively why the 
hardness of the electrons is not a strong function of radius 
if dark matter annihilations create the haze. Ignoring the 
boundaries of the diffusive region, the number of electrons 
at a given point in our dark matter models is given by 



n(E,r) = 



2B{av) [, 3 , 



M 2 

DM 



J dV p(|r + r'|) 2 G(r'\X{E)), (14) 

where p(r) is the dark matter halo profile, and the second 
line approximates the halo profile as the r" 1 inner cusp of an 
NFW profile with scale radius r s . The last line evaluates this 




Figure 5. Same as Figure \3\ but instead for four dark matter 
models. The red solid curve is a model with Mdm = 100 GeV, 
K = 10 28 cm 2 s _1 ,i= 0.3, and B = 10. The green dot-dashed 
curve is the same as the solid curve but with an inner power-law 
slope of —1.2, the blue short-dashed curve is the same as the solid 
curve but for a 10 times more massive WIMP and B = 800, and 
the magenta long-dashed curve is the same as the solid curve but 
with a 3 times larger value for Kg. 



expression in the limit A<r. When r > A, the hardness of 
the electron population or n(Ei,r) /n(i?2, r) essentially does 
not depend on r. This is in contrast to our other models, 
where the dependence of the hardness on r is exponential 
when r > A because n(E,r) oc G(r\X(E)). In the opposite 
regime where A > r, all haze models have a hardness that is 
a weak function of radius. The haze electrons cannot be in 
the latter regime at all radii and still reproduce the observed 
haze profile. 



4 COMPARISON OF MODELS 

Table [T] summarizes the parameters of 6 models for the 
haze emission that were motivated by the scenarios de- 
scribed in the previous section. The parameters of these 
models have been selected to yield a decent fit to the 
iDobler fc Finkbeinerl l|2008h measurement at 23 GHz in the 
inner 20°. Models A and B are two dark matter models for 
the haze (Section [33}, models C and D are models where it 
is due to a single event (Section 13. 4|) . and Model F is where 
the haze is due to a Galactic Center population of super- 
novae (Section 13. 3p . In Model E, an unknown Galactic Cen- 
ter population possessing a hard spectrum with 7i n t = 1.5 
is responsible for the haze (Section 13. 3. 3|l . 

Figure [5] compares the microwave intensity profile (left 
panels) and microwave hardness (right panels) of the 6 mod- 
els against that of the haze measurement. This comparison 
is the same as in Figure except for two differences. First, 
the haze measurement in the 41 GHz channel is also com- 
pared with the models (bottom panels). The shaded region 
that brackets the estimate for the systematic errors is not 
included in the bottom left panel because these systemat- 
ics are large at 41 GHz. Instead, only the profile from the 
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Table 1. Parameters of the different haze models plotted in Figure[6] All models assume U T + C/b = 10 eV cm -3 and = \/2/3 10 £tG. 
Models C-F assume that the input cosmic ray spectrum has sharp cutoffs at 1 and 1000 GeV. The "line-style" column provides the color 
and line-style used in Figure [6] 



model 


description 


Tint 


E tot or E tot (cgs) 


K (cm" 2 s- 1 ) 


<5 


line-style 


A 


100 GeV WIMP, B = 10 






10 28 


0.3 


red solid 


B 


1 TeV WIMP, B = 800 






10 28 


0.3 


green long dashed 


C 


GC Event, 1 Myr ago 


2.5 


8 x 10 51 


2 x 10 28 


0.6 


black dot-dashed 


D 


GC Event, 1 Myr ago 


2 


5 x 10 51 


2 x 10 28 


0.6 


magenta dot-long dash 


E 


GC Population 


1.5 


10 38 


2 x 10 28 


0.6 


brown dotted 


F 


GC Supernovae 1 


2 


2 x 10 38 


2 x 10 29 


0.3 


blue short-long dashed 



The soft synchrotron has been subtracted from this model as described in Section 13.31 




8 (deg.) 6 (deg.) 

Fig ure 6. Comparison of haze models described in Table 1 with 
the iDobler fc Finkbeinerl l|200Sl ) measurement. Left panels: The 
solid black curve with dots is the measured intensity profile of the 
WMAP haze in the 23 and 41 GHz channels. The other curves 
are the intensity profile predicted by the different models. Right 
panels: The ratio of the intensity in the 23 and 33 GHz channels 
(top panel) and of the 23 and 41 GHz channels (bottom panel) 
for these models. The curves in which the spectrum softens with 
radius represent the continuous source models, those in which 
this ratio is not a strong function of angle represent the dark 
matter models, and those that harden with radius represent the 
burst models. The column of black dots represent the measured 
hardness of the different haze estimators as described in Figure 
3. For comparison, the leftmost column of colored dots at 11° 
are the ratio of the model intensities in these two bands averaged 
over 6 — 15° , where colors are chosen to match the colors of the 
corresponding model curves. Lastly, the cyan short-long dashed 
horizontal line is the intensity ratio expected for free-free emis- 
sion, and the dotted black horizontal line is the intensity ratio of 
the galactic soft synchrotron radiation where I 



fiducial CMB estimator in IDobler fc Finkbeinerl (|2008h is 
included. Second, there are colored dots in the right panels 
that represent the mean hardness for the model curves with 
the corresponding color. These colored dots are calculated 
by taking the ratio of the intensity averaged between 6° and 
15° in the two channels. 

All of the considered model classes make different pre- 
dictions for the haze hardness. The qualitative trends for the 
hardness with radius are robust to the choice of diffusion and 
other model parameters. These trends rely on just the fact 
that more energetic electrons both diffuse more quickly and 
radiate their energy faster than less energetic ones. Models 
E and F (the continuous source models) predict that the 
spectrum should soften considerably with increasing angle, 
and Model F has the largest diffusion coefficient of the con- 
sidered models. Models C and D (the single burst models) 
predict the opposite, that the spectrum should harden with 
radius. Finally, models A and B (the dark matter models) 
predict that the spectral hardness should be almost constant 
as a function of angle. 

IDobler fc Finkbeinerl ||2008l ) measured significant emis- 
sion between 20° and 35° at the 1 kjy sr _1 level. All of our 
previous models for the haze emission begin to underpredict 
the emission compared to the haze measurement at > 15°. 
There are clai ms that the haze at these l arge angles may not 
be significant (|Cumberbatch et al]|2009h . If the fast decline 
in the haze intensity between 5° and 10° reflects the diffu- 
sion scale as we have so far assumed, it is difficult to explain 
the full extension of the haze with the previously considered 
models. However, if the diffusion scale corresponds to the 
cutoff at ~ 35° and the inner part of the haze owes to varia- 
tion in the magnetic field, then models with larger values for 
the diffusion length may be able to explain this exterior haze 
emission. Previous literature on the haze has taken different 
stances on whether to include in t heir fits the ou t er reg ions 
of the haze e mission. For e xample , iHooper et al.l (|200S ) did 
not, whereas ICholis et all l|2009bl ) did. 

Figure [7] shows two simple models, a supernova and a 
dark matter model, in which the diffusion length has been in- 
creased to fit this exterior emission (by a factor of 2.2 for the 
supernova model and 5.5 for the dark matter model relative 
to the corresponding models in Table 1). The models in this 
figure use A = 2 x 10" 16 GeV s _1 and B± = ^/2/3 5 (j,G, 
which are smaller than our fiducial values, and they both 
assume 8 = 0.3 and Lhazc = 10 kpc. These models do not 
attempt to fit the inner part of the haze, but allowing for a 
radial dependence of B± would likely provide the freedom 
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20 30 
9 (deg.) 



Figure 7. Same as Figure [3] but for two models that have been 
tuned to fit the outskirts of the microwave haze emission. These 
models use A = 2 X 1CT 16 GeV s" 1 and B± = 5 /J,G, which 

are smaller than the fiducial values, and they both assume S = 0.3 
and Lhazc = 10 kpc. The thin blue dashed curves represent a 
supernova model with 7; nt = 2, Kq = 2 x 10 29 cm -2 s _1 , and 
E = 3 X 10 38 erg s —1 . The corresponding thick curves are the 
same, but subtracting off the soft-synchrotron component of this 
model as described in Section [33] The solid red curve represents 



a dark matter model with Kq 
100 GeV, and B = 8. 



6 X 10 28 cm" 2 



Mr 



for these models to accommodate the measured emission 
profile. To fit the outer regions of the haze does not re- 
quire significantly more energy than in our previous models 
(despite the smaller B-field): The supernova model assumes 
3 x 10 38 erg s" 1 (a 50% increase over Model F) and 
the dark matter model assumes B = 8 (comparable to what 
is taken in Model A) . The angular dependence of the spectral 
hardness between these two cases is not as different as the 
cases in which A is smaller, but the predicted mean value is 
still much different between the dark matter and supernova 
models. 

In all of our simplified models for the haze, we have 
made numerous simplifying assumptions. For example, we 
fixed the value of Lhaio- We find that changing Z/haio does 
not significantly affect our conclusions because Lhaio > A in 
all of our models. We have also assumed a single magnetic 
field value. The synchrotron specific intensity is proportional 
to the line-of-sight integral over no E 3 ' B 2 . Noting that at 
fixed frequency E oc B~ 1/2 , this implies I v ~ B (1+l)/2 . 
If the Galactic B-field has a scale he ight of 2 kpc, as i s 
commonly assumed to fit radio data (|Strong et al.l l2007l ) . 
this will act to suppress the emission at latitudes of b > 
2kpc/8kpc rad. 



5 INVERSE COMPTON EMISSION 

Inverse Compton scattering off of the hard electron popula- 
tion that is responsible for the haze emission can contr ibute 
to the galactic 7-ray emission fe.g.. lHooper et alJ : 2008V An 
electron up-scatters via the inverse Compton process a pho- 
ton of energy -E S eed and wavelength A scc d on average to an 
energy of 



E-, 



40 



1 /im 



E 



80 GeV 



GeV, 



(17) 



where this equation is evaluated in the limit r\ = 
E se edE I (m e c 2 ) 2 -C 1. For E and A SO ed of interest, this limit 
is not always strongly satisfied. 

Models for the interstellar radiation in the central kpc 
of the galaxy find that starlight at 1+ 2 7 fim contributes 
most of the energy density in radiation, while ~ 10% is 
at A » 100 um from th e repossessing of this light by dust 
ijMoskalenko et al1l2006l ). The majority of the emission that 
results in the 7-ray at _E 7 ~ 10 (50) GeV likely owes to 
the Compton scattering of starlight off of electrons with 
E ~ 40 (90) GeV. These are slightly higher energy elec- 
trons than those that contribute to the microwave emis- 
sion. Therefore, this 7-ray emission is sensitive to slightly 
higher energy electrons than the microwave synchrotron. 
This statement assumes that 7 > 2. As 7 — ► 2, a compa- 
rable contribution to the 7-ray spectrum results from the 
inverse Compton scattering of ~ 100 fim dust emission on 
10 times more energetic electrons. 

Interestingly, haze-like emission may have already been 
detected in the firs t year data from th e Fermi Satellite in 
10s of GeV 7-rays (|Dobler et al.ll2010h . This emission was 
extracted in a similar manner to how the microwave haze 
was detected: First, galactic templates for no emission, for 
inverse Compton from the soft synchrotron electron pop- 
ulation, and for a uniform extragalactic background were 
subtracted from the Fermi 7-ray sky map in different spec- 
tral bands. Surprisingly, this resulted in a residual emis- 
sion with similar sp atial structure to the haze. Specifically, 
iDobler et~a l. (2010) showed that this emission could be fit 
by a bivariate Gaussian with standard deviations of 15° 
in longitude and 25° in latitude and with central intensity 
[E-^e]-, w 5 x 10" 7 GeV cm" 2 s" 1 sr" 1 at 40 GeV. 



We now investigate what the Fermi haze in combination 
with the microwave data reveals about the electron popula- 
tion in the haze region. For a power-law electron distribution 
assuming 77 <C 1, these measurements constrain the spectral 
index of these electrons to be 



2.7 - 0.9 In 



[E-jIe]-, 



10" 6 GeVcm- 2 s- 1 sr- 



5 kjy sr 

23GH 



E~< 



A s , 



O.'.Un [ ^1 +1.3 In 

(7* / \ 30 GeV 1 fim 



1 + 0.4 In 



30 GeV l^trn 



(18) 



where we have assumed for simplicity B± = \/2/3 10 ^iG 
and [/* is the energy density in starlight, which we have 
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assumed is concentrated at the wavelength of A scc dEf] How- 
ever, we find that a more realistic spectral distribution has 
a minor effect on the inferred value of 7. The logarithms 
evaluate to zero in equation (|f 8[) if values are taken for the 
Fermi and WMAP hazes that are consistent with measure- 
ments, yielding a soft spectrum with 7 ~ 2.7. Equation (| ISp 
primarily illustrates that there are significant astrophysical 
uncertainties inherent in the predictions of the 7-ray bright- 
ness of the the microwave haze electrons. For example, a 
factor of 2 uncertainty in Ub or [_E 7 /e] 7 results in an uncer- 
tainty of « 0.6 in 7. 

In addition, the inverse Compton spectrum alone can 
also be used to probe the spectrum of the ele ctron popula- 
tion. The measurement bv lDobler et al.l (|2010T ) of the Fermi 
haze emission finds roughly EjIe ~ E~ ' 2 at 40 GeV, with 
hints of a cutoff at higher energies. Such a slope only re- 
quires 7 ~ 3.4 (for 77 <C 1) and is not evidence for a hard 
electron spectrum. 

Different models for the microwave haze would pro- 
duce different trends in the spectral hardness of their 7- 
ray emission. If a single event produced the haze, it must 
have occurred At < 4 x 10 5 [10 eV cm~ 3 /((7 r + Ub)] X 
[40 GeV/-B 7 ] 1//2 yr ago for energetic enough electrons to be 
present to emit at E y for A scc( j = 1 /im. In addition, the 7- 
ray haze from an event would have a spectrum that hardens 
with radius, in contrast to the other haze models. If instead 
the haze owes to a Galactic Center population (such as su- 
pernovae), the inverse Compton spectrum will be harder to- 
wards the Galactic Center again because the most energetic 
electrons would have diffused the shortest distances. 

For dark matter models, the hardness in the 7-ray emis- 
sion should be fairly independent of radius from the Galac- 
tic Center because hard electrons are being produced in situ 
and do not have to diffuse from the Galactic Center. Most 
WIMP dark matter models also produce tens of GeV 7- 
rays more directly as annihilation byproducts. These direct 
byproducts should have an emission region that is less spa- 
tially extended than the inverse Compton emission from the 
haze electrons owing to diffusion effects. 



6 CONCLUSIONS 

We have investigated several models for the WMAP haze 
emission. We argued that free-free emission is unlikely to be 
the source of this radiation, in agreement with previous stud- 
ies. Therefore, it is most likely that the haze is synchrotron. 
If it is synchrotron, the haze indicates that the spectrum of 
electrons between a few and 20 GeV is harder towards the 
Galactic Center. 

We showed that it is difficult to produce such a hard- 
ened source of Galactic Center electrons with the standard 
disk population of supernovae. In addition, we argued that 
a disk population of pulsars (which has been cast by the lit- 
erature as the astrophysical alternative to the dark matter 

13 Equation 1 1181) is derived using the relation 

[EI e ]r = g| Ub N(E r ) 

[EIe]g E%U*N(E a )' 
where subscript R is for the radio haze, G is for the 7-ray, and 
N(E) oc E 1 is the electron spectrum. 



annihilation model) is also unlikely to produce a hardened 
spectrum just towards the Galactic Center. The bulk of this 
paper investigated different models that may achieve suffi- 
cient hardening. We showed that a Galactic Center popula- 
tion of supernovae and a diffusion-hardened electron distri- 
bution may be able to create the haze emission. We found 
that ~ 1000 Galactic Center supernovae per Myr could be 
consistent with present constraints on the energetics, profile, 
and (marginally) the spectral hardness of the haze if a rel- 
atively large value for the diffusion length is assumed. This 
supernovae rate is consistent with predictions for its value 
in the central couple hundred parsecs. 

A transient event in the Galactic Center that releases 
~ 5 x 10 51 erg in cosmic ray electrons in the past few Myr is 
another possibility that can reproduce the properties of the 
haze, even with a relatively soft spectrum (7i nt ~ 2 — 2.5). It 
is unlikely that a single stellar event can source such ener- 
gies, and, thus, such an event would be almost certainly as- 
sociated with Sagittarius A*. Intriguingly, there is evidence 
for enhanced star formation in the Galactic Center « 2 Myr 
ago. This activity could have been accompanied by a brief 
episode of accretion onto the SMBH, potentially producing 
an accelerated population of electrons. 

Synchrotron radiation from the byproducts of dark mat- 
ter annihilati ons is the most exciting explanation f or the 
WMAP haze l|Finkbeinerll2004bl ; iHooper et alj|2009h . This 
explanation is naturally able to account for the measured 
microwave hardness of the WMAP haze emission and its 
radial extent. This is true even if the diffusion length of 
~ 30 GeV electrons is much smaller than the radius of the 
haze, in contrast to the other haze models. A relatively light 
particle (< 100 GeV) is required in order not to invoke large 
boost factors in the annihilation rate beyond theoretically 
expected values. However, even for a 100 GeV WIMP, we 
require a boost factor of ~ 10 unless the inner halo density 
profile is steeper than a power-law with index —1.2. 

Electrons with energies that radiate in the microwave 
and 7-ray are expected to lose significant energy as they dif- 
fuse a distance comparable to the extent of the haze. Because 
of how these losses occur in the different models, all of the 
considered models produce differing trends in the hardness 
of the microwave and 7-ray spectra as a function of angle 
from the Galactic Center. These trends do not rely on the 
specific details of galactic cosmic ray propagation model. In- 
stead, all of these trends owe to the theoretically expected 
and empirically confirmed fact that more energetic parti- 
cles diffuse more quickly and cool faster than less energetic 
ones. Previous observations in the microwave and 7-ray are 
not yet able to constrain radial trends in the haze hardness, 
but future observations are projected to constrain them with 
high precision. 

We illustrated these trends usin g simple models that 
were t uned to be consistent with the iDobler fc Finkbeinerl 
(2008) haze measurement. In particular, dark matter annihi- 
lation models for the haze emission differ from astrophysical 
models in that a lot of the energy is input even at several 
kpc from the Galactic Center. Because of this difference, the 
hardness of the synchrotron and inverse Compton emission 
in this model is not as strong a function of galactic latitude 
as in the astrophysical models. This observable may allow 
one to distinguish this model from the others. 

Future observations have the potential to significantly 
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constrain potential explanations for the WMAP haze. In 
particular: 

• Radio maps at frequencies higher than 408 MHz would 
probe electrons that are more akin to the population that 
emits in the WMAP bands, and such maps could test the two 
electron population model for the microwave synchrotron 
emission. In fact, such data already exists at 1420 MHz 
l|Reich et alj|200lh . 

• Planck temperature data (and possibly more sophis- 
ticated techniques applied to the WMAP data) can pro- 
vide a better measurement of the haze profile as well 
as a measurement of the hardness as a function of an- 
gle. Planck's improved measurement relative to WMAP 
stems from the added channels at higher frequencies, which 
break degeneracies between the haze and other foregrounds 
i|Dobler fc Finkbeiner|[2008l ). 

• Future CMB polarization measurements can search for 
a polarized haze. While a recent analysis of WMAP po- 
larization data d id not detect a hard electron component 
l|Gold et al.ll2011T ). the haze should be present at some level 
if it owes to a hard population of synchrotron electrons. 

• 7-ray observations with Fermi are sensitive to slightly 
more energetic electrons than those that emit in the mi- 
crowave. Future 7-ray measurements have the potential to 
constrain the hardness of these electrons as a function of 
angle and thereby haze models. 

We thank Greg Dobler for providing the data for the 
iDobler fc Finkbeinerl (|2008l ) measurement of the haze, and 
Lars Hernquist, Philip Hopkins, Eliot Quataert, Avi Loeb, 
Ryan O'Leary, and Jennifer Yeh for useful discussions. We 
would especially like to thank Greg Dobler and David 
Spergel for useful comments on the manuscript. MM is sup- 
ported by the NASA Einstein fellowship. 
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